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Abstract. Abundances for a— elements and Fe in about 150 field subdwarfs and early subgiants with accurate 
parallaxes and kinematic data are used to discuss the run of abundance ratios in metal-poor stars in the solar 
neighborhood. Based on kinematics, we separated stars into two populations: the first one has a positive velocity 
of rotation around the galactic center, and it is likely to be related to the dissipational collapse of the galaxy; the 
second one has either negligible or negative rotational velocity, and it is likely related to an accretion component. 
The two populations show a large overlap in metallicity. However, they show distinct chemical properties. For the 
first population we found that there are close correlations (with small scatters around) of the rotational velocity 
with metallicity and with the Fe/a abundance ratio: this might be a signature of a not very fast collapse of the 
progenitor clouds, with enough time for a significant contribution by SNe la, although this result needs to be 
confirmed by a 3-D/non-LTE study. On the other side, the second population exhibits a larger scatter in both 
the above mentioned relations, and on average, a larger Fe/a ratio at a given metallicity, suggesting a larger 
scatter in ages. We argue that the lack of stars with moderate rotational velocities and high Fe/a abundance 
ratios is due to the short merging time for protogalactic clouds with prograde motion, while the presence of a 
group of counter-rotating stars with this characteristics indicates a much longer typical lifetimes for protogalactic 
fragments having such a motion. Finally, we found that perigalactic distances correlate with the Fe/a abundance 
ratios better than the apogalactic distances. 
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1. Introduction the surrounding interstellar medium is still rather limited 

' (Kay et al. 2002). Local Group galaxies, and in particu- 

Galaxy formation theory must account for the wide di- , r, n , , ., , 

. J . lar our own Galaxy, represent a unique opportunity to 

vcrsity oi observed structures, from dwarf irregulars to , , , . , 

J ' ° tune these models. Recently, given the improved capa- 

grand design spirals, from low brightness dwarfs to hum- r , • , i • , ln , 

° n i , i bilities ot new generation computers, wnicn allow us to 

nous ellipticals. This involves at least two critical ingre- , , ,. au a: a 

, . . carry out very high particle resolution runs, the hrst nu- 

dicnts: the formation of dark halos and global formation . , . , ,. £ n n , , r 

° merical simulations ol disk dominated galaxy formation 

of stars. The hierarchical clustering is currently the most . <■ nn . , , , . , r , i , 

° . m a fully consistent cosmological framework have been 
successful paradigm of structures formation. In the past c , , a i j p n v. j ™ n o al r a i nnan 

, , ..... . . . performed (bamland & Gerhard, 2003; Abadi et al. 2002; 

two decades numerical simulations have become an lm- „ , , , r, nn0 N m, , 

. Governato et al. 2002). lhere are, however, many unre- 
portant tool for exploring the detailed predictions of these , , , , a , , , , n , , 

r , , . . , ^ , ,. solved problems affecting both dark matter and galaxies m 

models (Weinberg et al. 2002). However, extending these , £ , , a , ■ onnl c .„ p D onnl i 

v & ' ' to such a framework (btemmez, 2001; bilk & Bouwens, 2001). 

studies to understand the formation of luminous compo- tw-ijj-j- r.i.v.ni u i j aj- 

, \ Detailed studies of the Galaxy can help our understanding 

nents of galaxies is difficult since our knowledge of the „ , , , . , , . , ,. ,. 

° . . . ° . . oi the complex processes involved m baryons dissipation 

star formation process and its interaction icedback) with / , ^ „ 01 , u , u of , no i 

1 v ' (see e.g. the review by freeman & Bland-Hawthorn, 2002). 
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atively uniform, isolated protogalactic cloud shortly after 
it decoupled from the universal expansion. In this scenario, 
gradients of metallicity and age with respect to kinemat- 
ics, as represented by rotational velocity or eccentricity, 
should be present (see e.g. Larson 1974, 1976), and the 
morphology of a galaxy should be related to the activity 
of star formation in early epochs compared to the collapse 
time. Angular momentum conservation causes the evolu- 
tion of the residual gas from an original roughly spherical 
distribution toward a disk. This picture was challenged 
by Searle & Zinn (1978). They found that the galactic 
globular clusters have a wide range of metal abundances 
independent of their radius from the galactic center. Their 
findings pointed out that individual protogalactic clouds 
might have had a significant independent chemical and dy- 
namical evolution before or during their merging into the 
main galactic body. Relics of this independent evolution 
may be traced by careful and extensive studies of chemi- 
cal and dynamical properties of the earlier galactic stellar 
generations. The ideas of galaxy formation via hierarchi- 
cal aggregation/merging of smaller units from the early 
Universe well agrees with the Searle and Zinn scenario of 
the galactic halo formation. However within this general 
framework, the role of the Milky Way is still to be ade- 
quately clarified. The Milky Way is a spiral galaxy located 
in a region of relatively low overdensity, at the edges of the 
Local Supercluster, with a quite significant bulge, a rather 
extended thin disk, a thick disk (Gilmore et al. 1989) and 
a halo. The presence of a dominant stellar (and gaseous) 
disk suggests that dissipational collapse, the main mecha- 
nism in the formation of spiral systems, had played a very 
important role. Moreover there are several observational 
facts that do not fit well in a monolithic dissipational col- 
lapse scheme. First, stars in the disk do not show the ex- 
pected run of increasing metallicity with increasing age; 
also there is a lack of metal-poor stars (the so-called G- 
dwarf problem). Both these unexpected features can be 
understood by assuming that there is a continuous infall 
of metal-poor material onto the disk. The origin of this 
material is not clear (although a number of interstellar 
clouds at high latitudes have been observed since a long 
time, starting with IRAS satellite). Second, the real dif- 
ferences between the different components of our Galaxy 
are still to be properly understood. Gratton et al. (1996, 
2000) and Fuhrmann (1998) have shown that the ratio be- 
tween a— elements (almost exclusively produced by short 
living core collapse Supernovae- SNe hereinafter) and Fe (a 
rather large fraction of which is produced on much longer 
time scale by thermonuclear SNe) indicates a clear dis- 
tinction between the epoch of formation of the thin and 
the thick disk, with a phase of low star formation in be- 
tween. This has led several authors to construct two-infall 
models, that reproduce quite well these observed features 
(see e.g. Chiappini et al. 1997). However, the reason for 
this hiatus in the star formation is still unclear. It may be 
due to either a merger event in the early history of our 
Galaxy, perhaps related to the formation of the bulge of 
our Galaxy, or the result of a galactic wind caused by a 



large number of SNe explosions in the later phases of the 
thick disk formation (or perhaps, to both these mecha- 
nisms). 

On the other side, the thick disk seems to have proper- 
ties that distinguish it from the halo too: in particular, it 
seems much more chemically homogeneous (see Nissen & 
Schuster 1997). It should be noted that the halo itself is of- 
ten considered to be composed of two distinct populations, 
one with some rotation, and the other either not-rotating 
at all, or even counter-rotating (Norris 1994). The second 
of these populations probably traces accretion events that 
occurred later in the galactic history; we are even actually 
observing one of these merging fragment, the Sagittarius 
dwarf Galaxy (Ibata et al. 1995). The relation between 
the rotating component of the halo and the thick disk is 
unclear, it may well actually be a single population with 
a continuum of properties. 

In this paper we present a discussion of the abundances 
of a— elements and Fe for a group of 150 metal-poor stars 
with accurate Hipparcos parallaxes, whose analysis was 
presented in a previous paper (Gratton et al. 2003, here- 
inafter Paper I) . High precision abundances and kinemat- 
ics were obtained for these stars: this allows us to discuss 
some properties of the the halo. In particular, we do not 
find any discontinuity between the properties of the rotat- 
ing component of the halo and the thick disk. They agree 
fairly well with predictions for a dissipational collapse of a 
large fraction of the primordial galaxy. On the other side, 
the not rotating component of the halo has clearly distinct 
chemical properties, in agreement with scenarios where it 
results from accretion of protogalactic fragments with in- 
dependent star formation history and chemical evolution. 

The paper is organized as follows: in Sect. 2 we recall 
the most relevant points of Paper I (sample selection and 
derivation of kincmatical data and chemical composition) ; 
in Sect. 3 we present the abundances, while in Sect. 4 we 
correlate these results with kinematical data and recall 
the most important biases present in our sample. Finally, 
in Sect. 5 we discuss our results in the framework of the 
early evolution of our Galaxy and present our conclusions. 

2. Sample and analysis 

Details about the analysis are described in Paper I. 
Briefly, we considered a sample of 146 mostly metal- 
poor stars with accurate parallaxes tt (At: /tt < 0.2) from 
the Hipparcos catalogue (Perryman et al. 1997); accurate 
space velocities could be derived for these stars, combin- 
ing Hipparcos proper motions and distances with radial 
velocities from the literature. To keep the analysis as ho- 
mogeneous as possible, we did not consider stars evolved 
off the base of the subgiant branch. High precision equiv- 
alent widths (with errors smaller than 3 mA) from high 
S/N, high resolution spectra were available for these stars, 
either from our own observations (using UVES at VLT, 
SARG at TNG, or the 2.7m McDonald telescope), or from 
the literature (Nissen & Schuster 1997; Fulbright 2000; 
Prochaska et al. 2000). 



Gratton R.G. ct al.: Abundances in metal poor-stars 



3 



The stars were subdivided into three groups according 
to kinematic criteria based on the results of appropriate 
Galactic orbit calculations. The groups considered are: 

1. A rotating inner population, with a galactic rotation 
velocity larger than 40 km s _1 and R max < 15 kpc, 
where i? max is the apogalactic distance. This popula- 
tion includes part of what usually is called the Halo 
as well as the Thick Disk, and can be substantially 
identified with the dissipative collapse population of 
Norris (1994), and with the population first identified 
by Eggen et al. (1962). Hereinafter, we will refer to this 
population as Dissipative Collapse Component. The 
reason for our choice is that we are not be able to find 
in our data any clear discontinuity between the prop- 
erties of the rotating part of the Halo and those of the 
Thick Disk. 

2. A second population composed of non-rotating and 
counter-rotating stars. This population includes the 
remaining part of what is usually called the Halo, 
and can be substantially identified with the population 
due to accretion processes first proposed by Searle & 
Zinn (1978). As we will discuss in the next sections, 
this population (as a group) has a chemical composi- 
tion clearly distinct from that of the first one, likely 
due to a different origin. We will call these Accretion 
Component Stars. 

3. Finally the Thin Disk, which also has clearly distinct 
chemical composition from the first population, as first 
showed by Gratton et al. (1996, 2000), and Fuhrmann 
(1998). The discontinuity in chemical composition be- 
tween the Dissipative Component and the Thin Disk 
is likely due to a phase of low star formation that oc- 
curred during the early evolution of our Galaxy (see 
also Chiappini et al. 1997). For such a component: 

v/Z2 mx + 4 e 2<0.35, (1) 

where Z max is the maximum height above the galac- 
tic plane (in kpc) and e the eccentricity of the galac- 
tic orbit of the star. With respect to the definitions 
of Paper I, we dropped the lower limit to [Fe/H] 
for the thin disk population. In this way we have 
a pure kinematical definition of the various popula- 
tions. Two stars were attributed to the thin disk rather 
than to the thick disk by dropping this criterion: HD 
116316 ([Fe/H]=-0.71, [a/Fe]=0.19) and HD134169 
(([Fe/H]=-0.84, [a/Fe]=0.28). While the impact of this 
change on the discussion is modest, we emphasize that 
stellar populations are defined on the basis of statis- 
tical criteria, and membership of individual stars may 
be questioned. 

A part from model uncertainties (which mainly affect 
absolute values, but only marginally the relative ranking), 
the largest source of errors in our kinematic is the error in 
the transverse motion related to the parallax. The mean 
quadratic error in the parallaxes of the stars of our sample 
is of An/n — 0.095. Since the average transverse motion 



is of 141 km s , the typical error in the transverse mo- 
tion is of about 14 km s _1 . This is much larger than the 
typical error in radial velocities, that for the stars in our 
sample is a few km s _1 . In order to estimate errors in 
the kinematical parameters, we then set up the following 
Monte Carlo procedure. For each star we randomly ex- 
tracted a series of values for the observational parameters 
(parallaxes, proper motions and radial velocities) obtained 
from Gaussian distributions centered on the observed val- 
ues and with dispersions equal to the errors assumed for 
these quantities. Errors in the parallaxes are those from 
the Hipparcos catalogue; for the remaining quantities we 
assumed typical values of 1 mas/yr for the proper motions 
and 2 km/s for the radial velocities. With this new set of 
data, we run our orbit extraction parameters. We repeated 
this procedure 100 times for each star, and adopted as er- 
rors in the orbital parameters the standard deviation of 
the distribution of the values obtained from each extrac- 
tion. With this procedure we obtain that typical observa- 
tional errors on the kinematic parameters we consider in 
this paper are ±0.39 kpc in the perigalactic distance i? m i n , 
±0.028 dex in the logarithm of the apogalactic distance 
logi? max , ±0.046 in the eccentricity e, and ±12 km s _1 in 
the rotational velocity VJ-ot- Maximum errors (that is, er- 
rors for the star which is most sensitive to possible uncer- 
tainties in the observational parameters) are about twice 
as large. Such errors are small enough to have a negligible 
impact on our discussion. 

A quite standard abundance analysis was carried out 
using these data; however, we carefully kept the star-to- 
star errors as small as possible. Effective temperatures 
for the program stars were obtained using B — V and 
b — y colors (generally assuming that interstellar redden- 
ing is negligible); these were calibrated against effective 
temperature using the Ha profiles that were available 
for about a third of the stars in the sample. A compar- 
ison with temperatures derived using the Infrared Flux 
Method by Alonso et al. (1996) shows an excellent agree- 
ment, with only a few slightly deviating point, gener- 
ally due to stars in binary systems. The final adopted 
temperatures should have typical errors of about 50 K. 
Surface gravities were obtained from the location of the 
stars in the color-magnitude diagrams, and should have 
errors of about ±0.1 dex. Finally, microturbulcnt veloci- 
ties were obtained minimizing any trend of derived abun- 
dances from expected equivalent widths for Fe I, and have 
errors of about ±0.3 km s _1 . Accurate, updated line pa- 
rameters were considered: in particular, collisional damp- 
ing was considered using the recent treatment by Barklcm 
et al. (2000), that should be much superior to the classical 
method based on enhancement factors over the classical 
van der Waals broadening given by the Unsold formula. 
Moreover, oscillator strengths for Mg were carefully re- 
vised, and non-LTE corrections were taken into account 
for O and Na lines. On the other side, we do not find any 
evidence for significant departures from LTE for the Fe 
lines used in our analysis. 
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Fig. 2. Run of the abundance ratio of Fe to the average of 
the a— elements vs [a/H]. Upper panel: data for the Thin 
Disk stars; middle panel: data for Dissipative Component 
stars; lower panel: data for Accretion Component stars. 
The mean line through various metallicity bins for 
Dissipative component is overimposed in each plot. The 
typical error bar is shown in the lower left corner 

Finally, two features must be mentioned. First, our 
analysis uses standard plane-parallel atmosphere models 
by Kurucz (1994, with the overshooting option switched 
off). We found a number of small but significant discrep- 
ancies in our results that can be ascribed to limitations in 
these model atmospheres; such limitations can perhaps be 
overcome using 3-d atmosphere models, with a much more 
realistic treatment of convection, such as those recently 
used by Nissen et al. (2002). Second, a quite significant 
fraction of the program stars arc members of binary sys- 
tems. Analysis for such stars is less accurate than for single 
stars: in particular, significant errors can be present in the 
adopted effective temperatures. As evidenced in Nissen et 
al. (2002) and in Paper I, both these concerns are impor- 
tant in the analysis of some key elements, like Oxygen. 
However, the impact on the elements considered in this 
paper is much smaller: for instance, we did not find any 
significant difference between the [Fe/a] ratios for bona 
fide single stars and for stars known to be members of 
binary systems (at variance from what we found for the 
[O/Fe] ratio). 

3. The Fe/a ratio, and kinematics 

In Paper I we considered abundances of several a- 
elements. In the present discussion we will only consider 
four of these elements: Mg, Si, Ca, and Ti. The most abun- 
dant among a— elements, O, is not included because its 



Fig. 3. Run of the abundance ratio of the average of 
the a— elements to Fe vs [Fe/H]. Upper panel: data for 
the Thin Disk stars; middle panel: data for Dissipative 
Component stars; lower panel: data for the Accretion 
Component stars. The mean line through various metallic- 
ity bins for Dissipative component is overimposed in each 
plot. The typical error bar is shown in the lower left corner 

analysis is more difficult for a variety of reasons (only 
the IR triplet at 7771-74 A could be used for the pro- 
gram stars; derived abundances are very sensitive to er- 
rors in temperatures; they arc affected by significant de- 
partures from LTE; and finally they arc quite sensitive to 
the adopted model atmosphere, with significant changes in 
case plane-parallel atmospheres are replaces by 3-d ones). 
We simply averaged the abundance ratios to Fe given by 
the different individual elements (Ti abundances being ac- 
tually the average of abundances given by neutral and 
singly ionized lines). Fig. 1 shows the correlation exist- 
ing between overabundances (with respect to Fe) of the 
various a— elements, using Ca as reference element: we se- 
lected Ca for this exercise because the [Ca/Fe] is the abun- 
dance ratio most accurately determined from our data. 
Different symbols are used for stars belonging to different 
populations. Error bars for abundances of individual ele- 
ments were determined in Paper I (see Table 9 of that pa- 
per). For those elements that will be considered through- 
out this paper, they are ±0.05 dex for [Fe/H] and [a/H]; 
0.025 dex for the [Fe/a] ratio (hereinafter a represents the 
average of Mg, Si, Ca, and Ti abundances); and 0.069 dex 
for the [O/Fe], 0.055 dex for the [Mg/Fe], 0.048 dex for 
the [Si/Fe], 0.033 dex for the [Ca/Fe], and 0.030 dex for 
the [Ti/Fe] abundance ratios. 

A 45 degree degree line is a good match to the ob- 
served data in all panels of Fig. 1, save for O. In this case, 
the slope is clearly larger, as expected since virtually no 
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Fig. 1. Comparison between the overabundances with respect to Fe for the individual a— elements with those obtained 
for Ca. Filled squares are data for Dissipative Component stars; open square those for the Accretion Component stars; 
crosses are data for Thin Disk stars. Error bars are shown in the upper left corner of each panel 



O is expected to be produced by type la SNe. The best 
fit line based on 66 stars (all but CD— 80 328, that ap- 
pears to be discrepant) is [O/Fe] = (1.72 ± 0.18)[Ca/Fc] + 
(0.089 ±0.1 18), where the error in the intercept represents 
the scatter of individual points. However, as mentioned in 
Paper I, O abundances are more uncertain for binary or 
reddened stars. Once these are eliminated, the best fit line 
is [O/Fe] = (1.80 ± 0.23) [Ca/Fe] + (0.046 ±0.100). Offsets 
are present for the remaining elements. On average, the 
abundance ratio to Ca are 0.132 ± 0.008 dex (r.m.s. of 
0.090 dex) for Mg, 0.002 ± 0.007 dex (r.m.s. of 0.086 dex) 
for Si, and -0.057 ± 0.006 dex (r.m.s. of 0.074 dex) for 
Ti, with no significant difference between stars belonging 
to the dissipative or accretion components. Some of these 
systematic differences might be due to offsets in the anal- 
ysis; however, the systematic decrease of the ratio with 
atomic weight corresponds well to the expected relative 
yields for different classes of SNe, in particular to a sys- 
tematic increase of the contribution by thermonuclear SNe 
with increasing atomic weight among the a— elements; or 
to a systematic increase of the yields for the lighter ele- 
ments by core-collapse SNe with increasing mass. While 



rather small, the r.m.s. scatter around the mean relations 
are larger than predicted from the errors listed in Table 9 
of Paper I: note that what is relevant here are the errors 
in the abundance ratios to Ca, that are expected to be 
0.078, 0.055, 0.051, and 0.043 dex respectively for O, Mg, 
Si, and Ti. The r.m.s. scatters remain higher even if a few 
(two to five) outliers are eliminated: in this case they are of 
O.lOOdex for [O/Ca], 0.085dcx for [Mg/Ca], and 0.067 dex 
for both [Si/Ca] and [Ca/Ti]. Taking into account the er- 
rors in the [Ca/Fe] ratio, there is still room for some real 
star-to-star scatter in these ratios, of about 0.063, 0.055, 
0.026 and 0.051 dex for the [O/Ca], [Mg/Ca], [Si/Ca] and 
[Ti/Ca] ratios. While it is well possible that these residual 
scatter are due to (underestimated) observational errors, 
we note that again, there is a trend with atomic number 
that may be suggestive of a real nucleosynthetic origin. 
However, the conclusions of this paper will not depend on 
these considerations. 

In the remaining part of this paper we will neglect 
the differences existing between individual a— elements 
(hereinafter not considering O), and consider their average 
abundance. Fig. 2 displays the run of the abundance ratio 
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Table 1. Average [Fe/a] values for the dissipative com- 
ponent in various metallicity bins 



fa /HI 


No. 


< fa /HI > 


fFe /al 

|i c/ aj 


rms 


bin (dex) 


Stars 


(dex) 


(dex) 


(dex) 


1.3 


4 


—1.715 


—0 415 ± 040 


0.079 


-1.3 - -1.0 


8 


-1.125 


-0.339 ±0.024 


0.067 


-1.0 - -0.8 


7 


-0.869 


-0.333 ± 0.020 


0.053 


-0.8 - -0.6 


12 


-0.688 


-0.346 ±0.019 


0.065 


-0.6 - -0.5 


15 


-0.544 


-0.278 ±0.016 


0.066 


-0.5 - -0.4 


16 


-0.443 


-0.286 ±0.020 


0.078 


-0.4 - -0.3 


8 


-0.342 


-0.238 ±0.030 


0.086 


-0.3 - -0.2 


10 


-0.253 


-0.261 ±0.023 


0.074 


-0.2 - 


5 


-0.116 


-0.226 ± 0.027 


0.061 



between Fe and the average of the a— elements with the 
average abundance of the a— elements. Stars belonging to 
the different components are shown in the three panels of 
this figure. Fig. 3 displays similar plots, this time for the 
abundance ratio between the average of the a— elements 
and Fe with respect to the Fe abundance. 

Tablc^lists the average abundances of the [Fe/a] ratio 
in different metallicity bins for the dissipative component. 
There is a clear trend for increasing [Fe/a] ratios with 
increasing metallicity: the linear correlation coefficient is 
r 2 = 0.257 over 85 stars, which is significant at a very 
high confidence level. The following best fit line seems to 
reproduce data very accurately: 

[Fe/a] = (0.113 ± 0.021)[a/H] - (0.229 ± 0.071), (2) 

where the error on the constant term represents the r.m.s. 
scatter along the line. This scatter is about three times 
the internal error (0.026 dex). 

Part of this scatter is probably still due to analysis 
errors; however, we may use it to set lower limits to the 
mass of clouds undergoing independent chemical evolution 
during the formation of this component of the Milky Way. 
Following the same approach described in Carretta et al. 
(2002), we will assume for this exercise that the scatter 
around an average [Fe/a] ratio is due to the stochastic 
effect of a sampling of the initial mass function by a dis- 
crete number of SN progenitors. Using the same precepts 
of Carretta et al., we then obtain that typically about 
> 100 SNe contributed to the enrichment of each cloud. 
The exact value of this lower limit depends on the uncer- 
tain dependence of yields on the stellar initial mass. On 
turn, the value of the Fe abundance may be used to con- 
strain the dilution factor of the ejecta of the SNe, and from 
the number of SNe and the yields we may then get an esti- 
mate for the mass of the whole cloud, that is of > 10 5 M@. 
This value is not much different from what Carretta et al. 
found for much more metal-poor stars, where the r.m.s. 
scatter but also the dilution factor are larger. As noticed 
by the referee of Carretta et al. (Bruce Carney) , this value 
for the mass of the clouds is of the same order of magni- 
tude of the Jeans mass at that epoch, and also of the same 
order of the mass of globular clusters. This may be per- 
haps considered then as a typical mass for clouds having 



an independent chemical evolution over a broad range of 
initial conditions in our Galaxy. 

We overplotted the line representing the average values 
of the [Fe/a] ratios for the dissipative component also on 
the two other panels of Fig. 2. Thin disk stars (shown in 
the upper panel) have a much larger average value of the 
[Fe/a] ratio at a given metallicity; the mean residual with 
respect to the best fit line of Eq. (2) is 0.176 ± 0.014 dex, 
with an r.m.s of only 0.059 dex, only slightly larger then 
the internal error. This large difference in the [Fe/a] ratios 
is simply confirmation of the results obtained by Gratton 
et al. (1996, 2000) and Fuhrmann (1998) for O and Mg. On 
the other side, stars belonging to the accretion component 
show a much larger scatter than that observed for stars of 
the dissipative component (the r.m.s. is of 0.108 dex over 
41 star). An F-test shows that the difference between the 
scatter around the mean line for the dissipative collapse 
and the accretion components is significant at 99.9% level 
of confidence, while a Kolmogorov-Smirnov test indicate 
that the two distributions are different at 95.0% level of 
confidence. Also the average [Fe/a] ratios of the two pop- 
ulations are different, in the sense that the [Fe/a] values of 
the accretion component at a given metallicity are larger 
by 0.043 ±0.017 dex; this is significant at 97.6% probabil- 
ity level (Student's t test). This result confirm an earlier 
finding (from a much more limited sample) by Nissen & 
Schuster (1997). 

Finally, before interpreting the observed slopes in the 
[Fe/a] ratio as real, we should warn the reader that we 
performed an LTE analysis using 1-d plane parallel model 
atmospheres. 3-D and non-LTE effects are expected to de- 
pend on [Fe/H], and may affect abundances of Fe and 
a— elements in different ways. Hence, a closer non-LTE/3- 
D study should be undertaken before the important con- 
clusion about type la SNe affecting the chemical evolution 
of the dissipative component can be considered definitive. 



4. Correlation with kinematics and biases 

Correlating our abundances with kinematics helps to 
understand several important features of our sample of 
metal-poor stars. In Fig. 4 we displayed the run of the 
metallicity (from Fe) with various kinematical parameters: 
Perigalactic distance i? m i n (upper left panel); Apogalactic 
distance i? max (lower left panel); Galactic orbit eccentric- 
ity e (upper right panel); Rotational velocity around the 
Galactic center v m t (lower right panel) . Different symbols 
are used for stars belonging to the various populations 
considered throughout this paper. Similar plots, but with 
(average) metallicity from the a— elements are shown in 
Fig. 5; while in Fig. 6 we plotted the kinematical quan- 
tities against the abundance ratio between Fe and the 
a— elements. 

Before opening a discussion on the results shown in 
these figures, we must warn the readers about the biases 
that affect our results. There are mainly three biases: 
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Our sample is drawn from the Hipparcos catalogue. 
This is complete only down to approximately V ~ 7; 
fainter stars (the vast majority of those included in 
the present sample) were obtained by combining lists 
proposed by several groups. These included a rather 
large number of high proper motion stars, so that there 
is a strong kinematical bias in our sample, favoring 
objects on highly eccentric orbits, with low rotational 
velocities around the galactic center, and with either 
large apogalactic or small perigalactic distances; 
Among the stars in the Hipparcos catalogue, we essen- 
tially selected metal-poor stars, that is stars known 
to have [Fe/H]< —0.8 (although a few stars more 
metal rich than this limit are actually present in the 
sample). The analysis of Carretta et al. (2000) shows 
that our sample is nearly complete insofar stars with 
[Fe/H]< —1.5 and accurate parallaxes are considered, 
but it is far from being complete for more metal-rich 
stars. However there is a clear bias against metal-rich 
stars; 



3. Finally our sample only includes stars with accurate 
parallaxes; this limits the volume sampled to a small 
region around the Sun, with a radius of the order of 
100 pc, the actual value depending also on the absolute 
magnitude of the stars. This biases our results against 
stars in the bulge, and also against objects of the outer 
halo on orbits of low eccentricities. 

It is very difficult to properly quantify these biases. 
However their impact on the results should be carefully 
considered in the following discussion. Our results cannot 
be used to draw correlations between metallicity and orbit 
parameters, a basic data that is much properly presented 
elsewhere (see e.g. Sandage & Fouts 1987, Norris 1994). 
Here, we will concentrate on a comparison between distri- 
bution of abundance ratios among different populations, 
under the assumption that they are less affected by the 
above mentioned selection effects. 

We note that in the case of the Dissipational Collapse 
Component the rotational velocity is clearly correlated 
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with the overall metallicity: the linear correlation coef- 
ficients are r 2 = 0.101 and 0.161 if [a/H] and [Fc/H] are 
used respectively (over 85 stars, that is 83 degrees of free- 
dom). Both these values are significant at more than 99% 
level of confidence. However, as mentioned above these 
correlations are at least in part artifacts of the selection 
biases present in our sample. We think that the (closer) 
correlation existing between the rotational velocity and 
the [Fe/a] ratio is much more interesting. In this case, the 
linear correlation coefficient is r 2 = 0.369, over the same 
number of stars. This correlation is significant at a very 
high level of confidence; as a matter of fact, rotational ve- 
locity can be used to predict the value of the [Fe/a] ratio 
with an accuracy better than that obtained using overall 
metallicity: the relation is: 

[Fe/a] = (0.00109 ± 0.00016)y rot - (0.454 ± 0.065). (3) 

The fact that the [Fe/a] ratio is better correlated with 
the rotational velocity than the overall metallicity is very 
interesting, because we expect that the effects of selection 



biases are in this case only indirect; that is, they should 
affect this correlation only through their impact on metal- 
licity and rotational velocity. 



5. Discussion and conclusion 

Our large sample of metal-poor stars with accurate abun- 
dance analysis and kinematics allows a discussion of some 
of the chemical properties of the components related to the 
dissipative collapse (rotating component) and to accretion 
(either counter- rotating or not-rotating). When examining 
the plots of Fig. 4-6, we should remind that a— elements 
are mainly produced by core collapse SNe. They then fol- 
low closely the star formation history of the Galaxy. On 
the other side, since a significant fraction of Fe is produced 
on much longer time scales by thermonuclear SNe, its pro- 
duction is somewhat delayed with respect to that of the 
a— elements. Hence we expect that a Fe/a ratio increasing 
with the time is a signature of a lower star formation rate. 
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The first feature we notice in our results is that 
even if both the dissipative collapse and accretion com- 
ponents show a large overlap in metallicity, they show 
distinct chemical properties. In fact, while in both the 
cases there is a significant deficiency of Fe with respect to 
the a— elements, at a given metallicity the Fe deficiency 
is more pronounced, and the scatter is smaller, among 
stars belonging to the dissipational collapse component, 
in agreement with earlier results by other authors (Nissen 
& Schuster 1997). As far as the dissipational component 
is concerned, the Fe deficiency outlined above is a signa- 
ture of a star formation rate higher than that of the ac- 
cretion component. Moreover the smaller scatter noticed 
above (Fig. 4-6) indicates that they formed from a more 
homogeneous interstellar medium that the accretion com- 
ponent, may be well mixed over regions typically as mas- 
sive as 10 5 M . Finally for the dissipational component 
we derive a close correlation between chemical composi- 
tion and kinematics, as expected for a continuous metal 
enrichment during the dissipative collapse phase of the 



Galaxy. However, while star formation occurred quite vig- 
orously, it still left time enough for a raise of the [Fe/a] 
ratio from a value of [Fe/a]~ —0.40 at [a/H]~ —1.5, to a 
value of [Fe/a]~ —0.23 at [a/H]~ 0, where the evolution 
of this component interrupts. This trend is very clear, but 
needs to be confirmed by a 3-D / non-LTE study. As found 
by other authors (Gratton et al. 1996, 2000; Fuhrmann 
1998; Reddy et al. 2003), this component is distinct from 
the thin disk, more rich in Fe, although shows an overlap 
in [a/H]. The raise of the [Fe/a] ratio during formation 
of the dissipative collapse component suggests that this 
lasted long enough to allow about 1/3 of the Fe produced 
by thermonuclear SNe to be deposited in the interstellar 
medium (again, with the caveat that this result might be 
due to 3-D/non-LTE effects not considered in the present 
analysis). While the time scale for the evolution of the 
progenitors of these SNe is still not well defined (see the 
review in Nomoto et al. 1999), this result suggests that 
the collapse phase was perhaps not very short. 



10 



Gratton R.G. et al.: Abundances in metal poor-stars 



On the other side, the accretion component exhibits 
a much larger scatter in the [Fe/a] ratio: this is essen- 
tially due to the presence of a significant population of 
Fe-rich stars, that made up nearly a third of the sample. 
We prefer the notation Fe-rich stars rather than the usual 
one, a— poor, because it is more suggestive of the origin 
of the anomaly: these stars are in fact quite metal-rich, so 
that an interpretation of their anomalies as due to random 
fluctuations of primordial SNe is not very likely (it would 
require that these stars formed in small clouds with a typ- 
ical mass of 3 10 4 Mq). Moreover within this framework 
we do not expect a significant difference between the av- 
erage [Fe/a] ratio of the accretion and dissipative collapse 
populations, at variance with what we find here. It seems 
more plausible that these stars formed in an environment 
with lower star formation rate and a chemical evolution 
quite independent of that of the Milky Way, similar to 
that found in the dwarf spheroidal galaxies, whose chem- 
ical composition has been recently studied by Shetrone et 
al. (2003) and Tolstoy et al. (2003), that halos tend to 
have higher Fe/a than stars in the dissipative component. 
Thus these stars should exhibit a large scatter in their 
ages. 

As discussed in Sect. 4, selection biases can affect pre- 
vious suggestions. Nevertheless we wish to call the atten- 
tion of the reader to the last panel in Fig. 6, which shows 
the trend of the [Fe/a] ratio with rotational velocity. As 
discussed above, the [Fe/a] ratio is most likely an alias of 
time. While there is a clear bias in this diagram in the 
sense that stars with low rotational velocity are under- 
represented, there is no obvious bias at a given rotational 
velocity. Thus it is interesting to compare the distribution 
of the values of the [Fe/a] ratio for stars belonging to the 
accretion component with that for stars with moderate 
rotational velocities (40 < V TO t < 150km s _1 ). While in 
the first group roughly half of the stars have [Fe/a]> —0.3 
(and so are probably significantly younger), stars as rich 
in Fe are very rare in the second group. The lack of slowly 
co- rotating stars with high Fe/a abundance ratios in the 
solar neighborhood is expected from dynamical models of 
disk formation. However a significant population of slowly 
counter-rotating Fe-rich stars points toward a long lifetime 
of their protogalactic clouds, in particular longer than that 
of the dissipational component. Moreover the possibility 
that accretion component stars belong to the bulge and 
that their kinematics properties have been changed by vi- 
olent or two-body relaxation processes, is not appealing 
in view of their low metal abundance ([Fe/H< —1). 

Finally, we wish to point out that perigalactic dis- 
tances correlate with the Fe/a abundance ratios much 
better than the apogalactic distances. At first sight, this 
result is quite surprising, since we would expect that the 
apogalactic distance is a better signature of the region 
of formation of the star and that accretion component 
stars should form mainly in the outer halo of the Galaxy. 
However, it should be noticed that the sample here consid- 
ered is drawn from the solar neighborhood, a region clearly 
distinct from the outer halo. In order to observe in the so- 



lar neighborhood stars formed in the outer halo, highly 
eccentric orbits are requested; if such an orbit passes very 
close to the galactic bulge, it becomes chaotic, loosing 
most of the information about its initial conditions (i.e. 
apogalactic distances) . 
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